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We have made a systematic investigation of charge transport in ID chains of Josephson junctions 
where the characteristic Josephson energy is much less than the single-island Cooper-pair charging 
energy, Ej <C Ecp- Such chains are deep in the insulating state, where superconducting phase 
coherence across the chain is absent, and a voltage threshold for conduction is observed at the 
lowest temperatures. We find that Cooper-pair tunneling in such chains is completely suppressed. 

Instead, charge transport is dominated by tunneling of single electrons, which is very sensitive to 
the presence of BCS quasiparticles on the superconducting islands of the chain. Consequently we 
observe a strong parity effect, where the threshold voltage vanishes sharply at a characteristic parity 
temperature T*, which is significantly lower than the the critical temperature, Tc- A measurable and 
thermally-activated zero-bias conductance appears above T*, with an activation energy equal to the 
superconducting gap, confirming the role of thermally-excited quasiparticles. Conduction below T* 
and above the voltage threshold occurs via injection of single electrons/holes into the Cooper-pair 
insulator, forming a non-equilibrium steady state with a significantly enhanced effective temperature. 

Our results explicitly show that single-electron transport dominates deep in the insulating state of 
Josephson-junction arrays. This conduction process has mostly been ignored in previous studies of 
both superconducting junction arrays and granular superconducting films below the superconductor- 
insulator quantum phase transition. 

PACS numbers: 85.25.-j, 73.23.-b, 74.25.-q, 74.78.-w, 74.81.Fa 


I. INTRODUCTION 

For several decades now, chains and arrays of low- 
capacitance Josephson junctions have attracted much 
attention as systems in which quantum phase transi¬ 
tions can be studied^i^, and as many-body platforms 
that could enable novel quantum phases^ and topo¬ 
logically protected states^i^. However, there have been 
many more theoretical proposals along these lines than 
experimental works, and experiments have mostly been 
confined to several intriguing avenues of research con¬ 
cerning the superconductor-insulator transition^, the dy¬ 
namics of quantum-phase slips2r2£, metrological current 
standardsii, conjectured solitonic phenomenai^rii, and 
the use of Josephson-junction chains as superinductorsi^ 

The insulating state of superconducting junction ar¬ 
rays is located below a superconductor-insulator (SI) 
quantum phase transition, and is synonymous with the 
destruction of superconducting phase coherence across 
the array, and localization of Cooper pairs. The Cooper- 
pair insulator occurs when the characteristic Cooper-pair 
charging energy significantly exceeds the Josephson cou¬ 
pling energy, Ecp Ej. The understanding and engi¬ 
neering of charge transport deep in the insulating state 
presents difficult problems due to the competition be¬ 
tween various modes of transport which include both 
Cooper-pair and single-electron tunneling processe o^^d'^ , 
both of which occur in the presence of significant disor¬ 
der. 

Nearly all studies of charge transport in insulating ar¬ 
rays, however, start from a point of view where it is as¬ 


sumed that the low energy excitations that play a role 
in conduction are Cooper pairs. This is based on the 
assumption that the temperature is sufficiently low to 
ignore the presence of unpaired electrons. A recent cal¬ 
culation of the dc conductivity for arrays deep in the in¬ 
sulating state, based upon single Cooper-pair excitations 
and including weak disorder, was put forth by Syzra- 
nov et a/3. Their model proposes Cooper pairs as the 
sole charge carriers, and describes transport in terms of 
variable range hopping between adjacent islands as a re¬ 
sult of Josephson tunneling. In another theoretical study 
by Fistul et a/3, a model for Cooper-pair transport in 
the insulating state of 1 and 2 D arrays of Josephson 
junctions is applied to the interpretation of experimental 
data from granular superconducting films. Similar lines 
of thinking have been taken in the analysis of Josephson- 
junction array experiments, also based on the assump¬ 
tion that charge transport far below the superconduct¬ 
ing transition temperature, is predominantly carried 
by Cooper pairs2S3. It is important to realize that in 
most experiments, the contribution of quasiparticles can¬ 
not be ignored. This contribution has proven important 
in studies of qubit3 and Cooper-pair transistor^, but 
has only recently been studied theoretically in JJ chains^ 

We have made a systematic investigation of charge 
transport for temperatures ranging from 10mK to IK 
in ID Josephson arrays that are characterized by large 
charging energies and high junction resistances. In this 
regime the arrays are deep in the insulating state. We 
find that Cooper-pair transport is completely suppressed 
and charge transport proceeds via single-electron in- 
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TABLE I. Device parameters 


device 

junction 

area 

(fimf 

island 

volume 

(pmf 

Rj 

Ecp 

Ej 

A 

0.015 

0.0029 

248 kD 

150 peV 

2.7 peV 

B 

0.015 

0.0029 

312 kD 

170 peV 

2.2 peV 

C 

0.003 

0.0016 

786 kD 

680 fieV 

0.86 peV 


jection into the Cooper-pair insulator. Furthermore, 
we observe a strong parity effecli^^, with a well-defined 
cross-over temperature T*, at which the voltage in¬ 
jection threshold decreases sharply. The parity effect 
has been studied previously in superconducting single¬ 
electron transistors^ and in Cooper-pair boxes^i^ 
where it appears as a temperature crossover from 2e to 
le periodicity in normalized gate charge. In contrast, our 
array measurements show the parity effect has a global 
effect on charge transport through the whole array. The 
presence of ~ 1 thermally-excited BCS quasiparticle per 
island in the array significantly enhances the tunneling 
rates of single electrons through the chain, and simul¬ 
taneously destroys the insulating state of the array, as 
the voltage threshold for single-electron injection is sup¬ 
pressed to zero. 


II. DEVICES AND MEASUREMENTS 

We have fabricated ID chains of Al-AlOx-Al Joseph- 
son junctions having a length of = 50 junctions using 
electron beam lithography, followed by thermal shadow- 
evaporation and in-situ oxidation of A1 films, which have 
a thickness of 30 nm. We have focused on three arrays 
with slightly different properties, having in common large 
junction resistances, Rj » Rq = (2e)^//i and large 
Cooper-pair charging energies, Eqp = {2e)'^ j^lCj ^ i5j, 
where Cj is the junction capacitance, as described in 
Table 1. The arrays exhibit a Coulomb blockade at low 
temperatures, both in the superconducting state as well 
as in the normal state, which is obtained by suppressing 
the superconducting gap in the films using an external 
parallel magnetic field. 

The samples were bounded to a circuit board, mounted 
in a microwave tight Cu sample enclosure, and secured 
to the mixing chamber of a BlueFors LD400 cryogen-free 
dilution refrigerator with a base temperature of 10 mK. 
Each DC line was filtered from high frequency radiation 
using three meters of ThermoCoax®, thermally anchored 
at each stage of the dilution refrigerator, and having a 
measured low-pass cutoff frequency of ~ 1 MHz. The 
lines were additionally filtered using chip LC components 
on the circuit board, and at room temperature using 
low-pass LC filters. Several measurements of supercon¬ 
ducting single-electron transistors were made using this 
setup that clearly showed 2e-periodic stability diagrams, 
indicating negligible quasi-particle poisoning due to the 
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FIG. 1. SEM image of the center region of device A. The 
center row of islands forms a Josephson jnnction chain with 
well-defined tunnel junctions. The angle-evaporation process 
also produces shadow islands, which are not galvanically con¬ 
nected to the junction chain. 


presence of non-equilibrium quasiparticles in the mea¬ 
surement leads. 

We have characterized the devices using current- 
voltage measurements with voltage biases ranging from 
5 to over 50 mV, and with a DC current resolution 
as low as '^0.8 femtoampere for currents up to several 
nanoamperes. In addition, we have used a parallel mag¬ 
netic field Hji to continuously suppress the supercon¬ 
ducting gap A(i?||). From our measurements we have 
found that the gap depends on the parallel magnetic field 
through A(i?||) ~ A(0)(1 —as expected^S, with 
A(0) = 210± 10/ieV, and = 0.59 ± 0.02T. 

It is important to distinguish between large-scale 
current-voltage characteristics (or IVCs) and small-scale 
IVCs. In the presence of a superconducting gap, ‘large- 
scale’ refers to voltage biases that span the onset of di¬ 
rect quasiparticle tunneling due to pair breaking across 
every junction in the chain. This occurs for eV ~ 2NA. 
‘Small-scale’ bias in the superconducting state refers to 
voltages V that are a substantially small fraction of 
2NA/e, in which case conduction is also referred to 
as sub-gap transport. As the superconducting gap is 
suppressed below approximately twice the characteristic 
single-electron charging energy, Eq = l2Cj = EcpjV 
the transition from small-scale to large scale occurs at 
V cs EcN/2e, above which the Coulomb blockade is 
lifted across each junction. Our main results are pri¬ 
marily concerned with sub-gap transport; however, the 
large-scale IVC measurements are used to experimen¬ 
tally characterize the energy scales and intrinsic disorder 
of the junction chains. 

IVC and differential conductance {dl/dV) data are 
shown for the three devices of Table 1 in Figure 2. 
The asymptotic normal state conductance at large bias 
voltage V determines the chain average of the normal 
state tunnel resistance Rj per junction. From this, one 
can extract the average Josephson coupling energy us¬ 
ing the Abegaokar-Baratoff relation Ej = ^A{Rq/R j). 
The Cooper-pair charging energy Ecp = {2e)'^l2Cj 
for each device is found experimentally by extrapolat- 
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FIG. 2. Large-scale current (upper plot) and conductance 
(lower plot) measurements for devices A, B and C in in the 
superconducting state for magnetic field (solid lines), 

and for device B in the normal state obtained with i?||=0.7 T 
(dotted lines). The inset in the upper plot shows the small- 
scale IVC at T = 25 mK for device B at femtoampere current 
resolution. Conductance {dl/dV) measurements in the lower 
plot are normalized to {NRj)~^ for each device. 


ing its normal state IVC data from the linear regime 
at large voltage bias to find the zero-current intercept 
Kffset = Ne/Cj. IVC's in the superconducting state in 
zero magnetic field show a steep onset of direct quasi¬ 
particle tunneling occurring at F ~ 2NAo/e, allowing 
one to directly estimate Aq. From the measurements in 
Figure 2, all devices are found to have nearly the same 
superconducting gap, Aq = 210 ± lO^eV. 

Devices A and B show a large ‘BCS peak’ that develops 
at y ~ 2NA/e. This peak arises from the overlap of di¬ 
vergent BCS quasiparticle density of states (DOS) of the 
superconducting islands of the chain. The peak is broad¬ 
ened primarily by offset charges that randomly shift the 
chemical potential of individual islands relative to their 
nearest neighbors. The result is a misalignment of the 
divergent BCS island DOS’s across the chain, which in 


the absence of charge disorder would be aligned when 
the voltage bias across individual junctions equals 2A/e. 
The broadening and concomitant reduction of this BCS 
peak therefore gives a measure of the relative disorder of 
the chains due to offset charges and fabrication inhomo¬ 
geneity. 

As clearly seen in Figure 2, the BCS peak for device C 
is significantly more broadened as compared to devices A 
and B. In addition, there appears to be a random struc¬ 
ture in the IVC data for this device around the onset for 
direct quasiparticle tunneling. This is expected as de¬ 
vice C has a significantly larger Eqp compared to device 
A and B, which makes it much more sensitive to ran¬ 
dom offset charges. In addition, device C was fabricated 
in a separate processing run, using different lithographic 
development parameters, and therefore may have more 
intrinsic disorder due to reduced film and junction qual¬ 
ity. 

Small-scale IVC data for device B is shown in the in¬ 
set of Figure 2. Note that the current scale in the inset of 
Figure 2 is 5 orders of magnitude lower than in the main 
plot. There is a clear voltage threshold V) ~ 1.3 mV for 
the onset of femtoampere currents. Threshold voltages 
can be distinquished quite clearly in the differential con¬ 
ductance at lower temperatures, as shown in Figure 3. 
The region \V\ < \Vt\ shows a current blockade and a 
zero-bias conductance Go = [dl/dV)v=o that is identi¬ 
cally zero, or at most, lower than our measurement reso¬ 
lution, G™'" = 10“^^ We experimentaly determine 
the voltage threshold V) as the absolute value of the volt¬ 
age bias at which the conductance rises above 10“^^ 

(see lower plot in Figure 3), which is a factor of 10 greater 
than the measurement resolution. As seen in Figure 4, 
there is a characteristic temperature T* at which V) drops 
sharply to zero, and above which a measurable zero-bias 
conductance is observed. 

A subset of the temperature dependent IVC data for 
device B is shown in the lower plot of Figure 3 for i3|| = 0 
and temperatures 25, 300, 350, 400 and 450 mK. T* for 
device B in zero field is found to be 270 mK (see Fig¬ 
ure 4). The data at 300 mK and above show a conduc¬ 
tance peak around V = 0 that grows with temperature, 
but starts out with a small mini-gap that appears to be 
a remnant of the blockade region below T*. The con¬ 
ductance peak arises from the overlap of the BCS DOS 
from island to island across the chain, and only becomes 
evident when there are thermally excited quasiparticles 
occupying these states. (For a large Josephson junction 
such thermally excited quasiparticles give rise to a loga¬ 
rithmic singularity at V = 0 in the IVC at finite T)4S 


III. PARITY EFFECT 

The dependence of the measured threshold voltage V) 
on temperature for device A is shown in the upper plot 
of Figure 4 for several values of the parallel magnetic 
held, i?||= 0, 0.25, 0.35 and 0.45 T. It is evident that V) 
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FIG. 3. Differential conductance dl/dV for device B in zero 
magnetic field. Upper: grayscale image of dl/dV on a 

logarithmic scale, showing a temperature dependent thresh¬ 
old that vanishes at a well-defined temperature, identified as 
the parity temperature T*. Lower: slices of dl/dV for tem¬ 
peratures 25, 300, 350, 400, and 450 mK. For temperatures 
below T*, the threshold voltage is clearly distinguished as a 
nearly two orders of magnitude increase in the conductance, 
as shown by the 25mK data. The voltage threshold Vt is 
experimentally determined to be the absolute value of the 
voltage bias at which the conductance rises above 10“^^ 
which is a factor of 10 greater than the measurement resolu¬ 
tion = 


vanishes sharply at a specific temperature that depends 
on i?||. We argue that this behavior is a consequence of 
the parity effect for small superconducting islands. The 
ground state free energy for an odd number of electrons 
is higher than that for even numbers by the amount F = 
A - /cBTlnTVeff, where fVeff(T) « VpiO)y/2TrkBTA{T) 
is the effective number of states at finite T arising from 
integration over the BCS quasiparticle DOS, V is the vol- 
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FIG. 4. Upper plot: temperature dependence of the voltage 
threshold Vt for device A for i3|| = 0, 0.25, 0.35 and 0.45 T. 
Lower plot: magnetic field dependence of the parity tempera¬ 
ture T* for devices A (triangles), B (circles) and C (squares). 


ume of the island and p(0) is the density of states for the 
normal metal at the Fermi energy2^. The free-energy 
difference F for a single island vanishes at a cross-over 
temperature k^^T* = A/In Afeff(r*). Using the island 
volumes given in Table I, the experimentally determined 
Aq, and taking p(0) = 1.45 x 10^^ m“^ J“^, for the 
density of states for aluminum, one can compute the 
theoretically expected parity temperature for isolated is¬ 
lands. For islands such as those in devices A and B, we 
calculate T* = 260mK, and for the islands of device C, 
T* = 277 mK. One notices that for our device parame¬ 
ters, k'B.T* Ri A/9, which is much less than Tq rsI.SK 
for aluminum. 

We experimentally determine the array parity temper¬ 
ature T* as the temperature at which Vt passes through 
the voltage threshold observed in the normal state at the 
lowest temperature, which is indicated in Figure 4 (upper 
plot) by the horizontal dotted line. For = 0 we find 
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T* = 256, 269 and 294 mK, respectively, for devices A, B 
and C. The dependence of T* on is shown in Figure 
4 for the three devices, and is observed to be linear with 
hence linear with A(i3||) as expected. Given that 
A(i?||) is found to be nearly the same for all devices, the 
results imply a different for each device. This result 
can only partially be understood within the single-island 
picture of the parity effect. The island volume for device 
C is nearly twice that for devices A and B, (see Table 1). 
Island volume affects T* logarithmically, and therefore 
would only lead to a variation in T* of 17mK. However, 
the observed difference in T* between device B and C is 
25mK, and between devices A and C, 38mK. Devices A 
and B, fabricated in the same evaporation and located 
in the same EBL field, show a difference in T* of 13 mK, 
even though they have nominally the same island volume. 

The data suggest that the inferred Nes for our devices 
does not follow from single-island considerations alone, 
but could be explained if the effective BCS quasiparticle 
DOS for the system of coupled islands is modified due 
to charge transport processes. One could expect Ngg to 
depend on factors such as the charging energy, tunnel 
resistance, and offset charge and other disorder in the 
array. As can be seen in the large-scale IVC data of 
Figure 2 (discussed above), device A with the lowest T* 
also exhibits the largest BCS peak in the tunneling DOS, 
followed by device B, which has the second lowest T*. 

Although the parity effect in superconducting single¬ 
electron transistors has been well known for the last two 
decades, it has been almost entirely neglected in stud¬ 
ies of Josephson-junction arrays. A sole theoretical pa¬ 
per by Feigel’man et al^ pointed out the implications 
of the parity effect on the experimental search for the 
charge-unbinding transition in 2D junction arrays. Dis¬ 
cussions of the parity effect in experimental works on 
both 2D junction arrays and 2D disordered supercon¬ 
ducting films, however, are conspicuously absent. Our 
results are particularly relevant for the latter, as the 
theoretical models describing disordered superconducting 
films below the superconductor-insulator transition are 
based on a junction-array picture for the mesoscopically- 
structured hlmsi^. 

In voltage-biased single-electron transistors, where 
there is no threshold for conduction, the parity effect 
is observed by change in the periodicy of the gate- 
dependent current at finite bias^^i^l. In a Cooper-pair 
box, the parity effect can be seen directly by measuring 
the box charge and observing a transition from le to 2e 
periodic Coulomb staircases^Si^. In the array measure¬ 
ments presented here, however, the parity effect has a 
global effect on charge transport across the whole chain of 
junctions. For a single island, the parity-dependent free 
energy difference goes to zero precisely as the thermal 
expectation for the number of BCS quasiparticles {Nqp) 
no longer depends on the charge parity of the island^^. 
In our junction chains, it appears that the presence of 
^1 thermally-excited BCS quasiparticle per island effec¬ 
tively destroys the insulating state, by sharply removing 


the voltage threshold for single-electron injection. In ad¬ 
dition, the tunneling rates for single electrons through 
the chain are significantly enhanced due to the presence 
of quasiparticles. The precise microscopic mechanism un¬ 
derlying this phenomena is currently being investigated. 
Some simulation results on parity effects in arrays can 
be found in Cole et al^. Recent experimental results on 
a hybrid normal-superconducting transistor illustrate en¬ 
hanced charge tunneling due to a non-equilibrium quasi¬ 
particle distribution!^ 


IV. THERMALLY-ACTIVATED 
CONDUCTANCE 

The zero-bias conductance Go = {dI/dV)v=o was 
measured under applied parallel magnetic fields for tem¬ 
peratures ranging from the parity temperature up to IK. 
Go in all devices is found to follow an Arrhenius law 
for thermal activation, Go(T) = Goo exp(—A^/I cbT), 
where Ea is the activation energy, as shown in the up¬ 
per plot of Figure 5. The zero-bias conductance in the 
normal state at 0.6 T (not shown) continues to decrease 
above T“^ = 10, indicating an electronic temperature 
lower than 100 mK— . 

As a function of applied B-field, the activation energy 
for devices A and B is linear with as shown in Figure 
5, and appears to be equal to A(i?||). The data can 
be accurately fit to Ea = A(i3||) = Ao(l — i3j^/i?^||), 
with Ao = 214 ± 3/reK, and Bc|| = 0.59 ± 0.02 T. We 
find that the experimentally determined value of Ea in 
zero applied field agrees with an independent estimate of 
Ao gained from the large scale IVC data, eVqp/2N = 
210 ± 10/reU, where is the voltage that marks the 
onset of direct quasiparticle tunneling that occurring for 
eVqp/N ~ 2A. An activation energy that equals the 
superconducting gap can be easily understood because 
an energy of 2A is required to break a Cooper-pair. Since 
two independent excitations are created, the exponent for 
thermal activation is A rather than 2A— 

For device C, Ea varies randomly and somewhat ir- 
reproducibly with i3||, taking values between 250 and 
350 /reV. We attribute this to significantly larger dis¬ 
order present in device C, as inferred from the large scale 
IV data, and which is also consistent with the much 
larger Ecp for device C. In contrast to this, however, 
the activation exponent for the conductance evaluated 
at V = 1.5mV for device C shows nearly identical behav¬ 
ior to that of Go for devices A and B. This voltage bias 
is just outside the observed Coulomb blockade region of 
device C in the normal state, which is also relevant for 
unpaired charge carriers when the islands of the chain are 
superconducting. We conclude that transport above the 
parity temperature is set by thermally-activated quasi¬ 
particles in chains where strong charge disorder does not 
dominate. 

Thermally-activated transport in ID SQUID-arrays 
was reported recently by Zimmer et al— . The use of 
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FIG. 5. Upper plot: log Go vs. 1/T for device B for i?||=0, 0.3 
and 0.4 T. The solid lines are fits to an Arrhenius law. Lower 
plot: activation energy Ea vs. for devices A (triangles) 
and B (circles). The dot-dashed line is a combined fit for both 
devices to A(B||) = Ao(l — B^/B^||) yielding zero-field gap 
Ao = 214 ± 3 /reVand | = 0.59 ± 0.02 T. 


a SQUID geometry permitted tuning Ej in situ using 
a perpendicular magnetic field. To account for their 
measurements, Zimmer et al. assume a zero-bias con¬ 
ductance that is the sum of two contributions: a flux- 
dependent part, as would be expected for Cooper-pair 
tunneling, and a flux-independent term that remains 
when Ej (Cooper-pair tunneling) becomes very small. 
With Ej suppressed to nearly zero, Zimmer et al. ob¬ 
serve an activation exponent of the order of the supercon¬ 
ducting gap. These authors interpret their measured Ea 
as a characteristic charging energy for localized Cooper 
pairs undergoing variable-range hopping, although they 
mention as an alternative explanation, thermally gener¬ 


ated quasiparticles. 

As far as we are aware, Zimmer et al. is the only 
reported measurement of thermally-activated zero-bias 
conduction in ID arrays in the superconducting state. 
Thermally-activated zero-bias conductance in 2D junc¬ 
tion arrays has been reported by two groups some time 
ago 35,36 _ These authors interpreted their results in terms 
of a so-called “core energy”, Ecore, which is the energy 
required to create an electron-hole pair on adjoining sites 
{e.g. by moving a single electron by one site), together 
with the induced polarization charge on neighboring is¬ 
lands. This model is known as the soliton model. In 
the superconducting state, one finds for a 2D system, 
Ecore = 2A -I- Ec/^^ where the first term come from 
breaking a Cooper pair, and the second is the electro¬ 
static energy for placing a single electron and hole on 
adjoining sites. Since two independent excitations are 
created, Ea = i?core/2 = A-|-£'c/4. While Tighe et alS^ 
found quantitative agreement with the core-energy model 
of localized dipoles, more detailed measurements by Dels- 
ing et al^ showed substantial deviations from this pic¬ 
ture. For large A/Ecp, Delsing et al^ interpreted 
their results as evidence for Cooper-pair/hole solitons, 
even though their measured activation energies showed 
a strong dependence on superconducting gap. Delsing et 
al^ also report thermally-activated conduction in a ID 
chain, but only in the normal state. 

For a ID chain with localized dipole excitations, one 
expects Ecore = 2A -I- Ec using the soliton model of 
Tighe et al. and Delsing et alM-^ , and therefore 
Ea = £lcore/2 = A -f Ecl‘2. = A -|- Ecp/S. Our results 
for ID chains show that Ea agrees more closely with A, 
with no additional term needed to account for the charg¬ 
ing energy of electron-hole pairs on adjacent islands. The 
localized dipole model ignores tunneling processes that 
effectively lower the core energy. As noted previously, 
above T*, a voltage threshold for conduction is no longer 
found. In summary, conductance above T* is consistent 
with the lack of an electrostatic threshold for both charge 
injection and activated transport. 


V. CONDUCTANCE BELOW THE PARITY 
TEMPERATURE 

Finally, we have measured conductance at 20 mK and 
above the threshold voltage V) as a function of the mag¬ 
netic field. Data for device A is shown in Figure 6 taken 
at a bias voltage, V = 4 mV. Here we have used the fit 
from Figure 5 to express B|| in terms of A. As the su¬ 
perconducting gap is suppressed by the magnetic field, 
the conduction Gj, in what we will call the “injection 
regime”, is clearly exponentially enhanced by the factor 
exp (—A/feeTefi). In contrast to the zero-bias conduc¬ 
tance, we find the effective temperature T^s = 340 mK, 
which is considerably larger than the zero-field parity 
temperature for this device, T* = 260 mK. This shows 
that charge transport below T* and above the voltage 
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FIG. 6. Magnetic field dependence of the conductance in the 
injection regime for device A at T = 20 mK. The conductance 
at V = 4 mV, Gi, is plotted against the experimentally deter¬ 
mined A(Bi^). The dot-dashed line is a fit to the expression 
InGi = InGj — A/fcBTefi, which yields = 340 mK. 


threshold occurs by injection of single electrons/holes 
into a non-equilibrium steady state, which shows a sig¬ 
nificantly elevated effective temperature. Future experi¬ 
ments are needed to address the detailed nature of this 
steady state and it’s relation to the voltage threshold for 
conduction, Vt, observed at low temperatures. 


VI. CONCLUSION 

In conclusion, we find that for ID Josephson-junction 
chains deep in the insulating regime, where Ej Eqp, 
there is a characteristic parity temperature T*, above 


which the insulating state is destroyed by thermally- 
excited BCS quasiparticles. Above T*, an observable 
zero-bias conductance appears and is thermally-activated 
with an activation energy equal to the superconducting 
gap. This can be understood most simply if charge car¬ 
riers are single electrons and holes rather than Cooper 
pairs. The situation is somewhat analogous to donor 
ionization in a doped semiconductor, although here the 
donors are the background of localized Cooper pairs, and 
the effective ionization energy is strongly renormalized 
due to the singular BCS density of states. Conduction 
at temperatures below T* occurring above the threshold 
voltage appears to be thermally-activated, with an ex¬ 
ponent equal to the ratio of the superconducting gap to 
an effective thermal energy, ksTeg. The effective tem¬ 
perature Teff is found to be significantly higher than the 
electronic temperature that would otherwise exist in the 
array. This indicates that a non-equilibrium steady state 
of unpaired charge carriers becomes established, enabling 
above-threshold charge transport below the parity tem¬ 
perature in the Cooper-pair insulator. Our results are 
also relevant to studies of disordered superconducting 
films, which are often modeled using a picture of weakly 
coupled superconducting islands. 
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